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Abstract 

We report on a search for the non-standard-model process u{c) -\- g ^ t using pp collision data 
collected by the CDF II detector corresponding to 2.2 fb~^. The candidate events are classified 
as signal-like or background-like by an artificial neural network. The observed discriminant dis- 
tribution yields no evidence for FCNC top-quark production, resulting in an upper limit on the 

production cross section (j(m(c) +g ^ t) < 1.8 pb at the 95% C.L. Using theoretical predictions we 
convert the cross-section limit to upper limits on FCNC branching ratios: B{t — > u-\-g) < 3.9 x 10~^ 
and B{t^c-\-g) < 5.7 x IQ-^. 

PACS numbers: 14.65.Ha, 13.85.Rm 
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In the standard model (SM) of particle physics the flavor quantum number of fermions can 
be changed by charged currents, i.e., weak interactions mediated by the exchange of a 
boson. Flavor-changing neutral-currents (FCNC) are absent at tree level, but do occur at 
higher order in perturbation theory through loop diagrams. These radiative corrections are 
further suppressed through the GIM mechanism In the bottom-quark sector the large 
top-quark mass alleviates the GIM suppression leading to FCNC decays with branching 
ratios at the level of 10~^, while in the top-quark sector FCNC decays are more strongly 
suppressed and occur only at the order of i3 ~ 10^^'^ to 10~^^ {2!, way beyond the current 
experimental sensitivity. Therefore, any evidence for FCNC in the top-quark sector will be 
a signal of physics beyond the SM. Enhanced FCNC effects can be realized in extensions of 
the SM, such as models with multiple Higgs doublets , supersymmetric models with 
R-parity violation [4], or top color-assisted technicolor theories In certain regions of 
parameter space of these models the branching ratio of FCNC decays can reach levels of 
10~^ to 10~^. But even with such an enhancement the detection of FCNC top-quark decays 
remains a very challenging task at the Tevatron: First, because one can only expect to 
reconstruct a few top quarks in these modes, and second, because the background for the 
most promising mode, t —>■ eg, is very difficult to discern from generic multijet production 
via quantum chromodynamics (QCD). It has therefore been suggested to search for FCNC 
couplings in top-quark production, rather than top-quark decay p, l7|. 

In this Letter we present a search for the non-SM single top-quark production processes 
m(c) + g ^ t. We do not consider a particular model, but perform a model-independent 
search based on an effective theory that contains additional flavor- changing operators in 
the Lagrangian 



9s^ u o^^U G% + g,^ c a^^^t G^^ + h.c. (1) 



A 2 ^ A 2 

Here Kfug and Ktcg are dimensionless parameters that relate the strength of the new, anoma- 
lous coupling to the strong coupling constant gs and A is the new physics scale, related 
to the mass cutoff above which the effective theory breaks down. The gluon fleld tensor 
is denoted GJ^^,, the A" are the Gell-Mann matrices, and cx^'^ = |[7'^,7'^] transforms as a 
tensor under the Lorentz group. The existence of FCNC operators allows the production of 
top quarks via u{c) + g ^ t, but also non-SM decays t —>■ u{c) + g. In the allowed region 
of parameter space for Ktug and Ktcg an experimentally favorable situation occurs. While 
the FCNC production cross-section of top quarks is in the range of several picobarns, the 
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branching ratio of FCNC decays is very small, and top quarks can thus be reconstructed 
in the SM decay mode t Wb. While u quarks are constituent quarks of the proton, c 
quarks, as needed for the process c + g ^ t, occur as sea quarks originating from a gluon 
splitting into a cc pair. In the SM, top quarks are either produced as ti pairs by the strong 
interaction or singly via the exchange of a virtual W boson. The pair-production process 
is firmly established experimentally with a cross section of about 7pb. Evidence for SM 
single top-quark production has been shown by CDF [8] and D0 [9i], yielding a cross section 
around 3 pb. 

Our analysis is the first one at the Tevatron searching for the 2-^1 processes u{c)+g t, 
while a previous D0 analysis [l^ has looked for 2 ^ 2 processes, such as qq — > tu, ug — > tg^ 
or gg tu, resulting in upper limits of Kfug/A < 0.037 TeV~^ and Hfcg/A < 0.15 TeV^^ 
at the 95% CL. FCNC couplings to the top quark involving the photon or Z boson have 
been constrained by the analysis of top-quark decays at the Tevatron [11], the search for 
e'^e~ — » tc/ tu reactions at LEP, see e.g. 121], and the search for ep ^ e + t + X reactions at 
HERA y, 



The analysis presented here uses pp collision data at y/s = 1.96 TeV collected by the 
CDF II detector |l^ at the Fermilab Tevatron between March 2002 and August 2007. The 
data set corresponds to an integrated luminosity of 2.2 fb^^. We select a set of candidate 
events in the t Wb iub topology based on the event selection used for the measurement 
of SM single top-quark production [8]. We require exactly one isolated electron with 
tra,.ve.=e ene.gy U Er > 20GeV o. one isolated muon with pr > 20 GeV/c. mining 
transverse energy > 25 GeV, and exactly one jet with pseudorapidity [16| \7]\ < 2.8 
and Et > 20 GeV. The jet is further required to contain a reconstructed secondary vertex 
consistent with the decay of a 6 hadron After all selection cuts we observe 2472 

candidate events. 

Background yields from diboson processes WW, WZ, ZZ, and ti production are pre- 



dicted using PYTHIA 
cross sections 



22 



2l| Monte Carlo samples, normalized to next-to-leading order (NLO) 



231 ]. SM single top-quark rates are estimated with simulated events from 



the tree-level matrix-element generator madevent 2J], subsequent showering with pythia, 
and normalization to NLO cross sections [251]. The processes with vector bosons {W or Z) 
plus jets are generated with ALPGEN 26|], with parton showering and underlying event simu- 



lated with PYTHIA. Using a compound model [8 



based on simulated events, theoretical cross 
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TABLE I: Predicted sample composition and observed number of jet events in 2.2 fb of 

CDF Run II data. 



Process 


Expected 


events 


Wbb, Wcc 


750.9 




225.3 


Wc 


622.3 




186.7 


Wqq 


769.9 




100.5 


tt 


12.3 




1.8 


QCD-multijet 


43.0 




17.2 


Diboson 


19.9 




2.0 


Z+jets 


26.6 




4.2 


SM single-top 


24.4 




3.6 


Total prediction 


2269.3 




434.3 



Observed 2472 



sections, and normalizations in background-dominated regions we predict the composition 
of the W+1 jet data set as given in Table [B Top-quark events produced via the processes 
u{c) + g t are simulated using the matrix-element generator toprex j^^], followed by 
parton showering with PYTHIA. For the event generation, the coupling constants have been 
chosen to yield a cross section of 1 pb, which correponds to the approximate sensitivity to 
the process with the data set we analyzed. By investigating kinematic distributions at par- 
ton level, we verified that the event kinematics do not depend on that choice of parameters 
within the range relevant for our analysis. Under the assumption that ntug = i^tcg the tug 
coupling contributes 0.94 pb and the teg coupling 0.06 pb. For a total FCNC top-quark cross 
section of 1 pb we expect a yield of 35.3 ± 5.3 events. 

For an efficient background rejection, we employ the same neural-network technology as 
used in the search for SM single top-quark production [s], [28|. Neural networks (NN) have 
the advantage that correlations between the discriminating input variables are identified 
and utilized to optimize the separation power between signal and background processes. 
The networks are developed using the neurobayes analysis package 29|], which combines 
a three-layer feed-forward neural network with a complex and robust preprocessing of the 
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input variables. The network infrastructure consists of one input node for each input variable 
plus one bias node, 15 hidden nodes, and one output node, which gives a continuous output 
in the interval [—1, 1]. We train the NN on the samples of simulated events listed above using 
a mixture of 50% signal events and 50% background events. The background composition 
is chosen in the proportions given in Table [H with SM single top-quark events included as 
background. In total, we use 14 variables that show significant discriminating power between 
signal and background. Variables derived directly from the four-vectors of reconstructed 
particles are the px and the 7] of the charged lepton, the px of the jet, the difference in azimuth 
angle between the jet and and between the lepton and as well as the Ai? between 
the charged lepton and the jet. The W-hoson candidate is reconstructed in its leptonic decay 
mode from the charged lepton and $rp applying the kinematical constraint M^^ = Myy = 
80.4GeV/c^. The two-fold ambiguity for the z-component of the neutrino momentum is 
resolved by choosing the smaller \pz,u\ solution. Based on the ly-boson reconstruction we 
define two input variables: the transverse mass Mt/u and rje^. We further reconstruct top- 
quark candidates by adding the jet to the reconstructed W boson and thereby define the 
following input variables: Mi^j, Mt/uj, the rapidity yi^j, and Qe-rji^j where Qe is the charge 
of the lepton. An additional input variable is the output of an advanced jet-flavor separating 
tool mainly developed to increase the sensitivity of the SM single top-quark searches 28|. 
To describe the event shape in general, we use the aplanarity of the reconstructed top-quark 
decay system [30 1. 

We apply the NN to the samples of simulated events and obtain template distributions 
of the network output for all physics processes considered. The template distributions of 
the most important background processes and the signal are shown in Fig. [JaJ As can 
be seen, the separation between FCNC top-quark events and SM single top-quark events is 
only marginal. The templates are weighted by their expected event yields and the resulting 
composite model is compared to the NN output distribution observed in collision data in 
Fig.Ifb) 



To measure the potential content of FCNC-produced top quarks in the observed data set, 
we perform a binned maximum likelihood fit of the NN output distribution. The effect of 
systematic uncertainties is parameterized in the likelihood function including the correlation 
of rate normalization effects and shape distortions of the template distributions. Uncertain- 
ties in the jet energy scale, 6-tagging efficiencies, lepton identification and trigger efficiencies. 
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the amount of initial and final state radiation, parton distribution functions, factorization 
and renormalization scale dependence, and Monte Carlo modeling have been explored and 
incorporated in this analysis. We integrate over all parameters describing systematic uncer- 
tainties in the likelihood function using Gaussian priors. The rate of Wbb and Wcc events is 
required to be positive, but otherwise unconstrained. Applying a prior probability density, 
that is zero if the FCNC cross section is negative and one elsewhere, we obtain the posterior 
probability density. No significant rate of top quarks produced by FCNC is observed and we 
set an upper limit on the cross section of 1.8 pb at the 95% C.L., which is in good agreement 
with the expected upper limit of 1.3 pb obtained from ensemble tests. The probability to 
obtain an upper limit higher than the observed 1.8 pb under the assumption that FCNC 
top-quark production does not exist is 28%. 



Using t 



effects 



31 



leoretical predictions of a{u{c) + g t), which include threshold resummation 



32l |. we convert the upper limit on the cross section into upper limits on the FCNC 



coupling constants at the 95% C.L. and find Ktug/^ < 0.018 TeV ^ assuming Ktcg = 0, and 



Ktcg/-^ < 0.069 TeV ^ assuming Kt^g = 0. Using predictions at NLO [33|, we also express 
these limits on the coupling constants in terms of limits on the FCNC branching ratios and 
obtain: 13{t ^ u + g) < 3.9 x 10"^ and B{t^c + g)< 5.7 x lO^^. 

For the first time we have explored the W+1 jet data set in search for top quarks produced 
by gluon-induced FCNC via the processes u{c) + g ^ t. No evidence for such processes is 
found, resulting in the most stringent limits on the branching fractions for FCNC top-quark 
decays. 
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FIG. 1: Distribution of the NN discriminant, (a) Discriminant shapes for the different physics 



processes normahzed to miit area. |(b)| The composite model is compared to the distribution 
observed in colhsion data. The inset shows the high NN-output region, where top-quark events 
contribute the most. 
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